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Theoretical  Determination  of  the  Heats  of  Formation  of 
Prospective  Strained-Ring  Rocket  Fuels 


J.A.  Boatz  and  J.D.  Mills 
OLAC  PL/RKS 

Propulsion.  Directorate 
Propulsion  Sciences  Division 
US-AFRL/Phillips  Laboratory 
Edwards  AFB,  CA  93524-7680 

ABSTRACT 

The  gas-phase  heats  of  formation  of  seven  unusual  strained  and  substituted  organic 
molecules  ranging  in  size  from  CgHg  to  C17H24N4O8  have  been  determined  using  the 
parallel  version  of  the  GAMESS  quantum  chemistry  code.  Molecular  energies  obtained 
by  a  number  of  methods,  ranging  from  the  semi-empirical  through  Hartree-Fock,  6- 
31(d)  geometry  optimization  and  force-constant  evaluation  with  MP2,  6-31(d)  single¬ 
point  energy  calculation,  have  been  combined  with  experimental  information  to  provide 
the  necessary  composite  thermodynamic  parameters.  In  addition,  a  variety  of  resource- 
tailored  strategies  for  employing  these  values,  some  relying  upon  investigation  of  the  target 
compound  alone  and  others  involving  suitable  isodesmic  reactions,  have  been  critically 
evaluated  and  compared  with  results  obtained  from  the  wholly  empirical  Benson  additivity 
rules.  Within  the  present  set  of  reference  molecules,  this  computation-free  approach 
compares  favorably  with  all  but  the  highest  level  theoretical  treatments. 


A.  Introduction 

As  a  measure  of  the  intrinsic  energy  content  of  a  substance,  the  standaxd 
enthalpy  of  formation,  the  energy  required  for  its  formation  from  the  elements  in 
their  standard  states  under  conventional  thermodynamic  conditions,  constitutes 
a  critical  parameter  for  the  evaluation  of  a  candidate  propellant’s  likely  perfor¬ 
mance.  Since  chemical  synthesis  of  even  small  quantities  of  novel  chemical  com¬ 
pounds  can  be  extremely  resource  intensive,  it  is  desirable  to  narrow  the  search 
for  new  materials  to  only  those  candidates  which  show  particular  promise.  In  or¬ 
der  to  be  of  use  to  propellant  chemists,  theoretical  predictions  of  molecular  heats 
of  formation  must  meet  the  often  conflicting  criteria  of  timeliness,  affordability, 
and  accuracy.  Thus,  when  encountering  a  set  of  perhaps  rather  speculative  pro¬ 
posed  fuels,  a  rapid  and  approximate  preliminary  screening  procedure  might  be 
productively  employed  in  advance  of  initial  experimental  investigations  and  more 
extensive,  and  presumably  more  accurate,  theoretical  treatments.  To  these  ends 
the  research  reported  herein  is  not  merely  motivated  by  the  quest  for  the  heats 
of  formation  of  a  few  specific  molecules  of  present  interest  but,  with  a  view  to¬ 
ward  future  needs,  also  gives  particular  emphasis  to  evaluating  and  comparing  a 
variety  of  established,  progressively  more  expensive  theoretical  methods  within  a 
diverse  set  of  strained  molecules  of  the  sort  of  interest  as  high-energy  propellants. 

Specifically,  this  work  reports  results  of  calculations  of  the  ideal-gas  en¬ 
thalpies  of  formation,  AHj  =  AHf®®  of  the  following  seven  compounds  made 
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or  proposed  by  synthetic  chemists  in  the  Propulsion  Sciences  Division  of  the 
Propulsion  Directorate,  Air  Force  Research  Laboratory. 


IX] 

1:  bicyclopropylidene 


4:  2,2,6,6-tetra(azidomethyl)- 
spiro[3.3]heptane 


X 

o 

$ 

o 

HO 

l:  1-cyclopropyl- 

V  V 

3:  2,6-dioxa- 

cyclopropan- l-ol 

spiro[3.3]heptane 

■N3  N02— v 

■n3  no2 — S 

V/vN — N02 

5:  2,2,6,6-tetra(nitromethyl)- 

spiro[3.3]heptane 


Nr^WW^N>  n°!^s/VW\^n°! 

N  — — N3  no, — 'Vn/vVn —  NOj 


6:  Higher  spiro  analogue  of  4 


7:  Higher  spiro  analogue  of  5 


In  addition,  the  specific  impulse,  a  parameter  derived  from  the  enthalpy  of  for¬ 
mation  and  more  directly  indicative  of  the  performance  of  a  molecule  in  a  rocket 
application,  is  reported  for  each. 

B.  Theoretical  Methods  for  the  Calculation  of  Heats  of  Formation 

There  follows  a  brief  survey  of  some  of  the  more  common  methods  devised 
for  the  prediction  of  molecular  heats  of  formation.  A  more  complete  account 
is  provided  in  a  subsequent  publication.1  These  approaches  may  be  broadly 
classified  in  a  number  of  ways;  for  present  purposes,  they  are  discriminated 
according  to  their  methodological  nature  or  scope  as  well  as  by  the  more  practical 
criteria  of  the  expense  of  the  calculations  which  supply  their  requisite  composite 
parameters. 

1.  Broad  Methodological  Classes 

Two  criteria  may  be  used  to  classify  the  methodology  of  a  determination  of 
the  molecular  heat  of  formation.  First  of  all,  some  approaches  involve  characteri¬ 
zation  of  only  the  molecule  of  primary  interest  (hereafter  the  “target  molecule”), 
in  contrast  to  those  to  be  discussed  subsequently,  which  seek  some  cancelation  of 
errors  by  utilizing  the  change  in  energy  of  a  reacting  set  of  molecules. 

Of  the  single-molecule  methods  considered  here,  the  second  criterion,  that 
of  additivity,  can  be  employed  to  further  discriminate  between  those  which  rely 
upon  the  approximate  transferability  of  energetic  contributions  of  chemical  moi¬ 
eties  such  as  atoms,  bonds,  or  functional  groups  (hereafter  labeled  as  ab  lateribus 
(literally  “from  the  bricks”))  and  those  which  conceive  of  molecules  as  complete, 
fully  interacted  systems,  indivisible  to  the  level  of  fundamental  quantum-chemical 

1  J.D.  Mills,  Technical  Report ,  in  preparation. 
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particles  ( “single- whole-molecule”  approaches).  In  these  latter  methods,  the  cal¬ 
culated  quantum-chemical  molecular  energy,  that  is,  the  energy  of  interaction 
of  the  composite  electrons  and  nuclei,  may  differ  from  the  enthalpy  of  forma¬ 
tion,  referenced  to  the  energies  of  the  elements,  by  several  orders  of  magnitude. 
Therefore  it  is  common  to  adjust  the  molecular  energies  by  a  variety  of  means  in 
order  to  correct  relatively  small  errors,  either  in  conjunction  with  or  in  advance 
of  incorporation  of  elemental  energies  and  thermodynamic  effects. 

An  alternative  means  of  correcting  minor  errors  in  whole-molecule  energies 
is  to  broaden  the  calculation  to  include  a  group  of  well-characterized  reacting 
compounds  containing,  in  some  sense,  similar  chemical  environments  or  groups 
(“reaction-based”  approaches)  and  to  take  advantage  thereby  of  a  partial  cance¬ 
lation  of  errors  between  the  calculated  energies  of  fully  interacted  reactants  and 
products.  These  approaches  then  separately  include  the  thermodynamic  effects 
of  non-zero  temperature. 

As  ordered  here,  the  methods  become  generally  more  complicated  and,  other 
things  being  equal,  are  expected  to  be  more  accurate.  They  also,  however, 
become  progressively  more  subjective,  as  for  example  in  the  selection  of  the  most 
appropriate  reaction  between  “most  similar”  reactants  and  products.  Each  of 
the  three  classes  already  described  are  more  definitively  characterized  by  the 
following  formulae. 

a.  Class  A -Ab  Lateribus  Methods 

The  enthalpy  of  formation  of  the  target  molecule  may  be  composed  as  a  sim¬ 
ple  sum  of  the  combined  energetic  and  thermodynamic  contributions  of  composite 
chemical  groups  as: 

AH*  =  Y,  AHflj  (1) 

j,  groups 

where  the  groups  may  be  operationally  defined  in  a  number  of  ways  so  as  to 
ensure  maximum  accuracy  and  transferability  between  different  molecules. 

b.  Class  B-Additively  Corrected  Single- Whole-Molecule  Methods 

Within  this  work,  methods  limited  to  the  fully  interacted  target  combine 
whole-molecule  energies  with  additivity-based  corrections,  as: 

AH}  =  E»„,  +  Y  AHj  (2) 

j,  groups 


somewhat  in  the  spirit  of  ab  lateribus  approaches.  These  undifferentiated  adjust¬ 
ments  are  designed  to  incorporate  thermodynamic  effects  as  well  as  ameliorate 
calculational  deficiencies. 

c.  Class  C-Reaction-Based  Methods 

In  the  present  research  isodesmic  or  “bond-preserving”  reactions,  those  in 
which  there  are  the  same  number  and  type  of  formal  bonds  in  the  sets  of  re¬ 
actants  and  products,2  are  used  to  determine  molecular  heats  of  formation  from 

2  W.J.  Hehre,  R.  Ditchfield,  L.  Radom,  and  J.A.  Pople,  J.  Am.  Chem.  Soc.,  92,  4796  (1970). 
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calculated  reaction  energies.  The  procedure  by  which  this  is  effected  is  per¬ 
haps  best  illustrated  with  a  specific  example.  The  third  target  compound,  2,6- 
dioxaspiro[3.3]heptane,  could  be  isodesmically  deconstructed  as: 

4CH4  +  oo  -4  2CH3-O-CH3  +  C(CH3)4,  (3) 

the  corresponding  reactions  for  the  remaining  target  molecules  being  suggested 
by  the  list  of  “non-target”  or  isodesmic-partner  compounds  listed  in  Table  (D). 
For  each  target  molecule,  the  ideal-gas  enthalpy  of  formation  is  then  constructed 
from:  - 

AHj  =  A(AH^)-  [AEeq+AEZp+A(H298K-HoK)j  (4) 

where  each  component  is  defined  below,  first  in  general  notation,  and  then,  as 
seems  illuminating,  in  application  to  the  example  reaction.  The  first  parameter: 

A(AH?)s  J>AHj(i') 

i;  molecules, 

1 

not  target 

=  2  AHj(CH3-0-CH3)  +  AHf(C(CH3)4)  -  4  AhJ(CH4)  (5) 

is  just  the  stoichiometrically  weighted  difference  of  known  heats  of  formation  of 
all  molecules  in  the  reaction  except  the  target.  The  quantity  enclosed  in  square 
brackets  in  equation  (4)  is  just  the  standard-temperature  reaction  enthalpy  and 
is  composed  of  three  parts: 

AEeq  ==  'y  Eeq(i) 

i,  molec. 

=  2Eeq(CH3-0-CH3)  +  Eeq(C(CH3)4)  -  4Eeq(CH4)  -Eeq(<<X» 
AEzp  =  'y  v{  EZp  (i) 

i,  molec. 

A(H298K_Hok)  =  (H298K(i)~HpK(i))  (6) 

i,  molec. 

which  are  similar  differences,  this  time  including  all  participants  in  the  reaction, 
of,  respectively,  the  energies  of  the  molecules  at  their  theoretical  equilibrium 
geometries,  shifts  due  to  zero-point  energy,  and  the  “heat  capacity  correction” 
which  incorporates  the  enthalpic  effects  of  the  standard  temperature.  To  the 
extent  that  the  structural  environments  are  similar  on  both  sides  of  the  reaction 
and  errors  in  a  calculation  are  similar  for  any  given  type  of  environment,  much 
of  any  error  is  expected  to  cancel  in  the  indicated  differences. 

2.  Calculational  Method 

Within  each  methodological  class,  specific  approaches  can  be  further  divided 
according  to  the  means  employed  to  supply  the  parameters  just  described.  In 
the  interests  of  determining  the  relative  value  of  each  method  in  a  resource- 
limited  environment,  it  seems  appropriate  to  organize  the  approaches  in  terms 
of  their  computational  expense.  Therefore,  labels  indicating  the  final  step  in  a 
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Table  (A) 

Calculation  Type  by  Computational  Expense 

Label 

Terminal  Calculation  Step 

I 

No  computation 

II 

Semi-Empirical  Geometry  Optimization  and  Hessian 

III 

Hartree-Fock  Geometry  Optimization 

IV 

Haxtree-Fock  Hessian  at  Minimum  Energy  Geometry 

V 

MP2  Energy  at  HF  Optimum  Geometry 

succession  of  calculations  (Table  (A))  will  be  used  hereafter,  in  addition  to  the 
methodological  class  index  already  described,  to  delineate  each  specific  approach; 
in  the  same  context,  the  label,  Exp,  will  be  used  to  denote  an  experimental  value. 

For  the  computations  performed  for  this  study,  the  semi-empirical  calcula¬ 
tions,  although  involving  the  same  effort  as  their  more  sophisticated  counterparts 
to  specify  the  geometry  and  prepare  the  input  files,  consumed  negligible  com¬ 
puter  time.  In  contrast,  each  succeeding  step  (III-V)  described  in  Table  (A)  was 
found  to  require  a  significant  and  roughly  comparable  additional  expenditure  of 
resources. 

3.  Details  of  Present  Methods 

a.  Method  IA-Benson  Group  Additivity 

Although  a  number  of  computation-free,  ab  lateribus  methods  have  been  de¬ 
vised,  this  work  considers  only  the  most  common,  the  Benson  group-additivity 
prescription.3  Based  wholly  upon  comparison  with  experiment,  representative 
enthalpic  contributions  of  chemical  groups,  each  generally  a  single  atom  discrim¬ 
inated  according  to  its  bonding  partners,  have  been  tabulated  and  may  be  quite 
readily  applied  to  many  molecules  with  classical  chemical  structures. 

b.  IIB  Methods-Semi-Empirical  Values 

A  number  of  computational  codes  directly  report  semi-empirical  enthalpies 
of  formation  after  shifting  the  fully  interacted  molecular  energy  by  atom-based 
adjustments  derived  from  experimental  and  theoretical  considerations.4  Results 
of  this  type  are  denoted  hereafter  with  the  label  IIB  followed  by  a  parenthetical 
indication  of  the  semi-empirical  Hamiltonian  used. 

c.  Method  IIIB-Ibrahim/Schleyer  Hartree-Fock  Atom  Equivalents 

Similarly,  Ibrahim  and  Schleyer,5  within  a  reference  set  of  molecules  for  which 
accurate  experimental  values  are  known,  found  an  optimum  set  of  atom-based 
shifts  which  yield  heats  of  formation  from  the  total,  Hartree-Fock  molecular 
energy.  These  factors  subsume  energetic  and  thermodynamic  effects,  as  well 
as  corrections  for  computational  deficiencies. 

3  S.W.  Benson,  Thermochemical  Kinetics,  Second  Ed.,  (Wiley,  New  York,  1976),  references  cited 
therein,  and  H.-D.  Beckhaus,  C.  Riichardt,  S.I.  Kozhushkov,  V.N.  Belov,  S.P.  Verevkin,  and 
A.  de  Meijere,  J.  Am.  Chem.  Soc.,  117,  11854  (1995). 

4  J.J.P.  Stewart,  J.  Comput.  Chem.,  10,  209  (1989). 

5  M.R.  Ibrahim  and  P.v.R.  Schleyer,  J.  Comput.  Chem.,  6,  157  (1985).  A  few  additional 
parameters  were  determined  by  the  present  author. 
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Table  (B) 

Origin  of  Parameters  (Eqns.  (5)  and  (6)) 
for  Reaction-Based  Methods  by  Calculation  Level 


Method 

AHf*(i') 

Eeq(i) 

Ezp(i) 

H298K(i)_H0K(i) 

lie 

Exp/G2 

SE 

SE(sc) 

Exp/SE(sc) 

me 

Exp/G2 

HF 

SE(sc) 

Exp/SE(sc) 

IVC 

Exp/G2 

HF 

HF(sc) 

Exp/HF(sc) 

VC 

Exp/G2 

MP2 

HF(sc) 

Exp/HF(sc) 

d.  Reaction-Based  Methods-Class  C 

The  sources  of  the  parameters  described  in  equations  (5)  and  (6)  for  each  of 
the  reaction-based  methods  are  summarized  in  Table  (B).  The  labels,  SE,  HF, 
MP2,  and  G2  denote  calculation  by  semi-empirical,  Hartree-Fock,  second-order 
Mqller-Plesset,  or  G2  methods,  respectively.  Semi-empirical  calculations  utilized 
the  MNDO,  AMI,  and  PM3  Hamiltonians,  as  will  be  indicated  by  a  parenthet¬ 
ical  addition  to  the  method  label,  and  Hartree-Fock  and  MP2  calculations  were 
performed  in  the  6-3lG(d)  basis.  The  label,  Exp,  followed  by  a  theoretical  level 
indicates  the  use  of  available  experimental  values,  either  direct  or  indirect,6  in 
preference  to  those  of  theory  and  (sc)  marks  theoretical  thermodynamic  values 
for  which  the  vibrational  frequencies  have  been  multiplicatively  scaled  (Hartree- 
Fock/6-31G(d)-0.897  and  semi-empirical:  MNDO-O.91,  AMI-0.95,  PM3-0.97,  as 
determined  by  the  present  author  by  comparison  of  the  zero-point  energies  of 
target  and  non-target  molecules  with  scaled  Hartree-Fock  values). 

4.  Computational  Considerations 

All  calculations  utilized  the  serial  and  parallel  versions  of  the  GAMESS  com¬ 
putational  code8  and  were  performed  on  IBM  RS/6000  workstations  at  AFRL  and 
tens  of  nodes  of  the  IBM  SP  machines  at  the  Maui  High  Performance  Computing 
Center  and  the  Aeronautical  Systems  Center,  Major  Shared  Resource  Center. 
The  largest  ab  initio  calculations  involved  molecules  with  no  symmetry  and  483 
basis  functions  and  consumed  several  node-months  on  the  IBM  SP. 

Quick  turn-around  and  high  through-put  are  enabled  by  the  extensive  paral¬ 
lelization  implemented  in  the  program.  For  example,  as  compared  with  sequential 
mode,  one  sample  MP2  calculation  ran  with  around  ninety-percent  efficiency  on 
eight  nodes. 

C.  Results  and  Discussion 
1.  Molecular  Geometries 

Hartree-Fock  optimization  indicates  that  molecules  1  and  3  have  the  re¬ 
spective  D2h  and  D2d  configurations  suggested  by  chemical  intuition  and  that 

6  Experimental  values  from  correlation  of,  J.B.  Pedley,  R.D.  Naylor,  S.P.  Kirby,  Thermochemical 
Data  of  Organic  Compounds ,  Second  Ed.,  (Chapman-Hall,  London,  1986),  pp.  89ff. 

7  J.A.  Pople,  H.B.  Schlegel,  R.  Krishnan,  D.J.  Defrees,  J.S.  Binkley,  M.J.  Frisch,  R.A.  Whiteside, 
R.F.  Hout,  and  W.J.  Hehre,  Int.  J.  Quantum  Chem.:  Quantum  Chem.  Symp.,  15,  269  (1981). 

8  M.W.  Schmidt,  K.K.  Baldridge,  J.A.  Boatz,  S.T.  Elbert,  M.S.  Gordon,  J.H.  Jensen,  S.  Koseki, 
N.  Matsunaga,  K.A.  Nguyen,  S.  Su,  T.L.  Windus,  M.  Dupuis,  and  J.A.  Montgomery,  Jr.,  J. 
Comput.  Chem.,  14,  1347  (1993). 
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these  minima  constitute  the  only  thermodynamically  relevant  geometries.  In 
contrast,  molecule  2  seems  to  have  no  symmetry  and  the  Hartree-Fock  poten¬ 
tial  energy  surface  exhibits  three  degenerate  pairs  of  closely  lying  (within  ~1 
kcal/mol)  enantiomeric  minima.  This  is  expected  to  modestly  degrade  the  qual¬ 
ity  of  the  thermodynamic  predictions  reported  herein,  as  they  assume  a  single 
harmonic  minimum. 

The  four  remaining  compounds  were  initially  subjected  to  an  optimization 
constrained  by  D2d  symmetry  with,  on  each  end  of  the  molecule,  the  nitrogen- 
containing  groups  directed  away  from  one  another  and  the  distal  methylene  hy¬ 
drogens  eclipsed  after  the  manner  of  propane.  Although  computation  of  the  full 
Hessian  revealed  only  real  frequencies  for  one  of  the  molecules,  the  remaining 
three  had  “imaginary  frequencies”  (magnitude  <  34  cm-1)  indicating  that  these 
are  not  equilibrium,  minimum-energy  structures.  By  sampling  a  number  of  ge¬ 
ometries  resulting  from  rotation  about  the  distal  C-C  and  C-N  bonds,  for  each 
molecule  an  almost  degenerate  pair  of  “constrained  minima”  of  D2  and  S4  sym¬ 
metry  were  found  to  have  an  energy  fully  4  to  7  kcal/mol  below  the  corresponding 
D2d  structures.  However,  vibrational  analysis  indicated  that  each  was  not  a  local 
minimum  on  the  full  potential  energy  surface.  Only  by  allowing  non-planarity  in 
the  spiro-linked  cyclobutane  rings  (as  in  cyclobutane  itself)  were  true  local  min¬ 
ima  discovered.  A  number  of  combinations  of  physically  justifiable  ring-bending 
distortions  starting  from  both  the  D2  and  S4  geometries  were  sampled  and  sev¬ 
eral  nearly  iso-energetic  minima  were  elucidated.  For  each  molecule  the  most 
stable  geometries  have  no  symmetry  and  lie  fully  5  to  11  kcal/mol  below  the 
D2d-constrained  optima.  The  principal  hope  that  these  indeed  constitute  global 
minima  rests  upon  the  progressive  diminution  of  the  stabilization  energy  with 
each  successive  type  of  distortion  from  D2d  and  the  extensive,  but  admittedly 
not  exhaustive,  sampling  of  the  alternatives.  In  any  case,  deficiencies  in  the 
thermodynamic  predictions  due  to  the  presence  of  multiple,  low-lying  minima 
and  an  anharmonic,  relatively  flat  potential  energy  surface  are  again  likely  to  be 
manifest. 

2.  Enthalpies  of  Formation 

Calculated  heats  of  formation  for  both  target  and  non-target  molecules  are 
presented  in  Tables  (C)-(E).  For  the  first  two  classes  of  methods,  those  limited  to 
a  single  molecule  (A  and  B),  results  for  the  target  molecules,  Table  (C),  may  be 
compared  with  those  for  the  non-target  molecules,  Table  (D).  As  evidenced  by 
the  root-mean-squared  deviations  from  either  experimental  values  (RMS(Exp)) 
or  (in  lieu  of  experimental  information)  the  highest  level  theoretical  predictions 
(RMS(VC)),  the  computation-free,  ab  lateribus  method  of  Benson  clearly  out¬ 
performs  the  other  tabulated  approaches.  This  appears  to  be  true,  not  only  for  the 
non-target  set  containing  relatively  small  and  simple  systems  expected  to  be  well- 
represented  among  the  reference  molecules  of  parameterized  empirical  methods, 
but  also  for  the  target  compounds  as  well.  Although  these  latter  systems  have 
well-characterized  classical  chemical  structures,  they,  not  unexpectedly,  manifest 
predicted  heats  of  formation  with  a  much  stronger  overall  variance,  either  by 
virtue  of  their  size  or  because  of  the  unique  character  of  their  strained  binding. 

The  predicted  target  heats  of  formation  displayed  in  Table  (E)  for  the 
reaction-based  methods  (class  C)  indicate  that  the  use  of  isodesmic  reactions 
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Table  (C) 

Target-Molecule  AH|  (kcal/mol) 

Ab  Lateribus  (A)  and  Single- Whole-Molecule  (B)  Methods 


Molecule 

IA 

IIB(MNDO) 

IIB(AMl) 

IIB(PM3) 

iiib 

1:1X1 

86.6 

61.8 

79.9 

74.0 

83.8 

2:  [>-|<] 

HO 

-8.5 

-12.4 

-0.4 

-5.4 

-9.1 

3: 

-26.9 

-62.0 

-36.8 

-44.6 

-30.2 

4:^OOC; 

320.6 

307.3 

342.7 

306.8 

281.1 

5- 

NOj«AAA-noj 

-19.8 

68.5 

1.6 

-26.1 

-27.4 

*ipooooc 

354.3 

313.8 

365.5 

316.2 

308.7 

DooooC 

13.9 

73.3 

24.9 

-18.1 

-1.4 

RMS  Dev.  vs.  VC,  RMS(VC) 

1  11.9 

_ 

52.9 

23.5 

13.2 

19.0 

N 

Ab  Lateribus  (A)  and  Sing 

Table  (D) 

on- Target-Molecule  AH|  (kcal/mol) 

le- Whole-Molecule  (B)  Methods  Compared  with  Experiment 

Molecule 

IAa 

IIB(MNDO) 

IIB(AMl) 

IIB(PM3) 

IIIBa 

Exp 

ch4 

-15.3b 

-11.9 

-8.8 

-13.0 

-17.0 

-17.8  ± 

0.1 

CH3N3 

66.4 

76.7 

69.9 

60.7 

71. 0C 

CH3N02 

3.3 

-9.9 

-15.9 

-17.9 

-17.8±  ( 

).l 

c2H6 

-20.4 

-19.7 

-17.4 

-18.1 

-20.1 

-20.0  ± 

0.1 

CH30CH3 

-43.6 

-51.2 

-53.2 

-48.3 

-43.6 

-44.0  ± 

0.1 

CH3CH2CH3 

-25.3 

-24.9 

-24.3 

-23.6 

-25.2 

-25.0  ± 

0.1 

cyclopropane 

12.8 

11.2 

17.8 

16.3 

13.6 

12.7  ± 

0.1 

C(CH3)4 

-40.3 

-24.6 

-32.8 

-35.8 

-40.2  ± 

0.2 

t-BuOH 

-75.1 

-64.3 

-71.6 

-71.3 

-68.8 

-74.7  ± 

0.2 

X 

-17.1 

-13.2 

-16.3 

-21.9 

-10.1 

-16.3  ± 

0.3 

RMS(Exp)  1  1.0  1  9-6  |  6.0  |  3.6  |  4.5  | _ | 

a  Absent  results  constitute  the  sole  value  in  a  parameter  reference  set  and  so  are  identical  to 
the  experiment  by  definition. 
b  From  bond  additivity  rules. 

c  G2  theoretical  result  from:  D.W.  Rogers  and  F.J.  McLafferty,  J.  Chem.  Phys.,  103,  8302 
(1995). 

tends  to  increase  the  quality  of  the  results  obtained  using  a  given  calculational 
level.  Further,  none  of  the  methods,  including  the  next  most  sophisticated  re¬ 
quiring  the  rather  expensive  Hartree-Fock  optimum  geometry  and  Hessian,  sig¬ 
nificantly  out-performs  the  Benson  ab  lateribus  method.  In  addition,  it  will  be 
noted  that  the  method  labeled  IIC(PM3)  seems  roughly  competitive  with  the 
Benson  rules.  Although  this  computationally  economical  semi-empirical  method 
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does  require  more  human  effort,  it  may  be  useful  in  cases  in  which  suitable  ab 
lateribus  parameters  are  not  available.  Its  utility  for  systems  lacking  classical 
chemical  structures  certainly  bears  further  examination. 

It  must  be  conceded,  however,  that  all  characterizations  of  target-molecule 
results,  heretofore,  rest  upon  the  assumption  that  the  most  expensive  method, 
VC  (or,  in  conventional  notation,  MP2/6-3lG(d)//HF/6-3lG(d)),  closely,  or  at 
least  best,  represents  experiment.  In  a  similar  study  of  a  variety  of  strained 
hydrocarbons  of  the  type  which  provide  the  framework  of  the  target  molecules, 
it  was  found  that  methods  essentially  identical  to  those  here  labeled  IVC  (or 
HF/6-31G(d)//HF/6-3lG(d))  and  VC  predicted  isodesmic  heats  of  formation 
with  RMS  deviations  from  experiment  of  10.4  and  2.1  kcal/mol,  respectively. 
These  values  seem  at  least  consistent  with  the  11.7  kcal/mol  deviation  of  the 
IVC  results  from  those  of  the  highest  level.  Further,  it  has  been  found10  in 
a  different  test  set  that  electron-correlation  treatments  of  a  sophistication  just 
beyond  the  MP2  can,  in  fact,  lead  to  reaction  energies  of  reduced  accuracy.  In  any 
case,  until  experimental  values  or  significant  additional  computational  resources 
become  available,  the  current  results  would  seem  to  provide  the  best  measure  of 
target-molecule  formation  enthalpies. 


Table  (E) 

Target-Molecule  AH| (kcal/mol)  by  Reaction-Based  Methods  (C) 

Molecule 

IIC(MNDO) 

IIC(AMl) 

HH 
1— 1 

a 

hd 

S 

CO 

IIIC(PM3)a 

IVC 

VC 

1 

61.4 

70.8 

82.1 

82.5 

75.1 

76.5 

2 

-19.7 

-9.0 

-12.2 

-9.9 

-9.9 

-10.2 

3 

-45.1 

3.3 

-25.4 

-27.4 

-29.8 

-28.8 

4 

318.9 

362.3 

332.5 

333. 9b 

329.3 

315.8 

5 

-22.4 

12.9 

-12.5 

-15.7b 

-18.8 

-34.7 

6 

312.0 

400.2 

348.0 

365. 0b 

358.0 

342.7 

7 

-31.0 

51.2 

1.5 

13. 7b 

8.4 

-8.6 

RMS(VC) 

17.7 

42.0 

11.6 

15.7 

11.7 

a  For  brevity  only  IIIC(PM3)  results  are  provided.  See,  J.D.  Mills,  Technical  Report ,  in 
preparation,  for  IIIC(MNDO)  and  IIIC(AMl)  values. 

b  Optimized,  semi-empirical  structure  differs  qualitatively  from  the  Hartree-Fock.  Thermody¬ 
namic  parameters  are  calculated  for  the  local-minimum  geometry  most  resembling  the  Hartree- 
Fock. 


3.  Specific  Impulse 

The  implications  of  the  theoretical  results  for  actual  propellant  performance 
may  be  illuminated  in  a  number  of  ways.  Most  simply,  a  candidate’s  mass- 
normalized  energy  content  seems  a  more  appropriate  measure  of  its  lifting  capa¬ 
bility  than  the  common  molar  value;  mass-adjusted  heats  of  formation  for  the 
highest  level  theoretical  results  are  given  in  Table  (F).  Further,  calculation11  of 

9  R.L.  Disch,  J.M.  Schulman,  M.L.  Sabio,  J.  Am.  Chem.  Soc.,  107,  1904  (1985). 

10  J.R.  Van  Wazer,  V.  Kello,  B.A.  Hess,  Jr.,  and  C.S.  Ewig,  J.  Phys.  Chem.,  94,  5694  (1990). 

11  C.  Selph,  R.  Hall,  C.  Beckman,  R.  Acree,  T.  Magee,  “Theoretical  ISP  Code,”  Computer 
Program,  Phillips  Lab.,  Edwards  AFB,  CA. 
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the  specific  impulse,  Isp,12  provides  a  more  direct  measure  of  the  probable  utility 
of  a  prospective  fuel.  Even  though  the  target  molecules  are  envisioned  as  solid- 
or  liquid-propellant  additives,  the  present  study  does  not  address  the  energetics 
or  properties  of  their  condensed  phases.  Therefore,  Table  (F)  is  limited  to  the 
specific  impulse  of  each  fuel  as  a  gas  combusted  with  liquid  oxygen  under  fixed 
reference  conditions.12  It  was  hoped  that  these  values  would  provide  at  least  a 
relative  measure  of  likely  performance;  under  the  same  conditions,  RP-1,  a  com¬ 
mon,  kerosene-based  hydrocarbon  fuel,  has  a  theoretical  specific  impulse  of  299.8 
sec. 


Table  (F) 

Target-Molecule,  Mass-Based  Enthalpies  of  Formation  (Method  VC) 
and  Specific  Impulse 

Molecule 

AH^cal/g) 

lap  (sec) 

f  1:  X 

955 

313.2 

2:  [>-|<] 

HO 

-104 

300.7 

3:  °00° 

-288 

295.5 

4:pOOC: 

998 

301.6 

-104 

287.8 

pooooc 

864 

301.8 

*  DOOOOC 

-21 

292.4 

Somewhat  surprisingly,  the  majority  of  the  proposed  fuels  seem  to  manifest 
relatively  poor  performance  or  only  modest  performance  gains  compared  to  RP-1. 
Additional  calculations14  reveal  that  their  specific  impulses  even  show  a  generally 
disappointing  relationship  with  those  of  the  non-target  isodesmic  counterparts- 
molecules  which  might  be  characterized  as  containing  the  same  functional  groups 
as,  but,  for  the  most  part,  lacking  the  energy-bestowing  strain  of,  the  targets.  It 
was  initially  supposed  that  the  relatively  unsaturated  target  molecules  might 
simply  lack  sufficient  hydrogen  to  produce  significant  quantities  of  favorable, 
light  exhaust  products  when  burned  with  oxygen  alone  and  that  they  might 
therefore  manifest  a  cooperative  enhancement  with  a  relatively  hydrogen-rich 
fuel.  However,  no  tripropellant  combination  of  liquid  oxygen  and  a  target  fuel 
with  either  liquid  hydrogen  or  RP-1  leads  to  any  predicted  performance  increase. 
(That  is,  in  each  case  the  optimum  specific  impulse  is  just  that  of  a  binary  mixture 
of  oxygen,  and  whichever  of  the  other  two  components  has  the  best  performance 
with  oxygen  alone.)  It  appears  that  condensation  energetics  and  densities  must 
be  determined,  or  at  least  sensibly  hypothesized,  in  advance  of  optimization  of 
realistic  specific  formulations  in  order  to  more  definitively  rank  these  candidates. 

12  G.P.  Sutton,  Rocket  Propulsion  Elements:  An  Introduction  to  the  Engineering  of  Rockets,  Sixth 
Ed.,  (Wiley,  New  York,  1992). 

13  Sea  level  expansion  with  1000  psi  chamber  pressure. 

14  J.D.  Mills,  Technical  Report,  in  preparation. 
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D.  Conclusions 


The  accuracies  of  a  series  of  progressively  more  expensive  and  sophisticated 
theoretical  methods  for  calculating  gas-phase  enthalpies  of  formation  have  been 
evaluated  in  a  specialized  reference  set  of  strained  and  substituted  molecules  in 
an  effort  to  justify  means  by  which  valuable  experimental  and  theoretical  re¬ 
sources  may  be  more  efficiently  focussed  on  candidates  of  particular  promise.  It 
appears  that  the  computation-free,  ab  lateribus  Benson  prescription  based  upon 
correlations  with  experimental  values  can  be  expected  to  predict  the  formation 
enthalpies  of  these  types  of  molecules  with  an  error  on  the  order  of  10  kcal/mol. 
As  an  alternative  to  very  lengthy  calculations,  facile,  semi-empirical  computation 
with  the  PM3  Hamiltonian  in  conjunction  with  isodesmic  reactions  may  provide 
roughly  comparable  results  for  the  relatively  rare  systems  to  which  the  simpler 
method  can  not  be  applied.  After  an  initial  screening  with  either  of  these  two 
approaches,  it  appears  that  the  progressively  sophisticated  theoretical  methods 
considered  here  do  not  yield  significant  improvements  in  accuracy,  at  least  until 
the  level  at  which  energies  from  a  Hartree-Fock-optimized  vibrational  analy¬ 
sis  and  MP2  single-point-energy  calculation  (MP2/6-31G(d)//HF/6-31G(d))  are 
combined  in  accord  with  an  appropriate  isodesmic  reaction. 

By  themselves,  theoretical,  gas-phase  heats  of  formation  constitute  parame¬ 
ters  determinative  of  a  simple  relative  specific-impulse  ranking  of  candidate  fuels 
under  reference  conditions.  Nonetheless,  condensed-phase  properties  and  more 
interaction  with  researchers  in  propellant  formulations  appear  to  be  necessary 
in  order  to  fully  actualize  the  potential  of  these  investigations  to  support  the 
effort  to  reduce  the  time  and  expense  required  to  develop  and  deploy  propellants 
capable  of  higher  performance  and  increased  payload. 
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